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Abstract 

We consider the time evolution of the density matrix p in a 2-dimensional com- 
plex Hilbert space. We allow for dissipation by adding to the von Neumann equation 
a term D[p\, which is of Lindblad type in order to assure complete positivity of the 
time evolution. We present five equivalent forms of D[p\. In particular, we connect 
the familiar dissipation matrix L with a geometric version of D[p\, where L consists 
of a positive sum of projectors onto planes in R^. We also study the minimal num- 
ber of Lindblad terms needed to describe the most general case of D[p\. All proofs 
are worked out comprehensively, as they present at the same time a practical pro- 
cedure how to determine explicitly the different forms of D[p\. Finally, we perform 
a general discussion of the asymptotic behaviour i — > oo of the density matrix and 
we relate the two types of asymptotic behaviour with our geometric version of D[p\. 
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1 Introduction 



Particle physics is not only a field where fundamental interactions are explored, but it has 
also become a testing ground for possible deviations from the quantum-mechanical time 
evolution. The time evolution of the density matrix p is given by the master equation 

^^=-^Hp + ^pH^-D[p], (1) 

where D[p] is a dissipative term which adds to the quantum-mechanical term on the 
right-hand side of Eq. ([l|). Such a term can emerge if the system under consideration is 
not fully closed but interacts weakly with the environment. In general the nature of such 
an interaction is unknown, but experimental data from suitable systems can be used to 
place bounds on the parameters of such hypothesized interactions. In the general case, 
the Weisskopf-Wigner approximation |l| allows to incorporate also unstable particles, by 
using non-hermitian Hamiltonians H = M — ir/2, where M and F > are hermitian. 

It is not only necessary that the time evolution (|I]) respects p{t) > Wt, but the 
stronger requirement of complete positivity [Q] seems to be a natural and physical concept 
(for the general structure of completely positive maps, see Refs. f^, Q). This concept 
is defined in the following way. Let us assume that the system under consideration is 
described by elements of the finite complex Hilbert space Ti (dim?-^ = d < oo) with 
time evolution p{t) = 7t(p) with p(0) = p. Considering in addition the finite-dimensional 
Hilbert space C^, one can extend the time evolution 7^ on to a time evolution 7„.j on 
H ® by defining ■jn-t = 7t ® 1. If the time evolution 74 derived from Eq. ([l|) with 
the dissipative term D[p] has the property that '~in;t{,p) > is valid for all times t, all 
n = 0, 1, 2, . . . and all density matrices on the space H. ® C", then 7^ is called completely 
positive. 

Complete positivity of 7f determines the general structure of D[p\ p, |[ (see also 
Refs. 0, 1, HITI). It has been shown by Lindblad § (see also Ref. Q) that 7* of Eq. (|I|) 
is completely positive if and only if D[p\ has the structure 

i=i 

= lt{[AhA,p] + [pA],A,]) , (2) 

where the operators Aj act on Ti. 
Using the relation 

Aj = Bj + Sj 1 with Tr Bj = 0, (3) 
the dissipative term (^) can be reformulated as 

D[p] = -t[AH, p] + D'[p] with AH='-j2 i^jB] - s]Bj) , (4) 

where D'[p] is obtained from D[p] by the replacement Aj Bj. This reformulation 
has the effect that part of D [p] is shifted into the quantum- mechanical term of the time 
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evolution (|T]), such that a new Hamiltonian H' = H+AH appears. Note that for hermitian 
operators Aj we have AH = 0. In the space of traceless operators on the Hilbert space 
H we can choose a basis {Fj \ j = 1, . . . ,cP — 1} with the property 



Tr (F/Ffc) = 5 



and expand the operators Bj as 



Then we obtain the expression 



'jk 



Bj = CkjFk . 

k=l 



(5) 



(6) 



D'[p] = E ( PB}] + [BjP, B]] ) = -- E ( [Fk, pF^ + [Fkp, F} 



kl=l 



(7) 



where {cki) is a positive matrix defined by 



Cki — E ^kjCi, 



The form (|^) of the dissipative term has been derived by Gorini, Kossakowski and Sudar- 
shan 0. It is equivalent to the Lindblad form (^). Thus the time evolution with H and 
D[p] is equivalent to the one with H' and D'[p]. 

In particle physics, searches for deviations from the quantum-mechanical time evolu- 
tion are going on in neutral meson-antimeson systems [K^K^ and B^B^) (for a list of 

P|, |2D|, |2T|) and neutrino physics (for a 



papers see, e.g. 



Refs. 
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The importance of complete positivity, 
in particular, in K^K^ and analogous systems, has been stressed in Refs. |16, jTS). For 
instance, only if 7^ is completely positive, the positivity of the time evolution 7j®7t in the 
space TC<^TC is guaranteed. This follows from the decomposition 'jt^lt = (7t® 1) (l®7t), 
where both factors are positive according to complete positivity. Forming the tensor prod- 
uct '-ft ® 7t is a method for implementing the 1-particle time evolution at the 2-particle 
level, which is an often used procedure, e.g., in the K^K^ system. For an example where 
7t is only positive but not completely positive, with ensuing non-positivity of 7^ ® '-ft, see 
Ref. m 



For simplicity we assume a hermitian Hamiltonian H from now on, but this assumption 
is irrelevant for our discussion of complete positivity; it only concerns the investigation 
of the asymptotic limit of the time evolution, where loss of probability due to a non- 
hermitian H as obtained in the Weisskopf-Wigner approximation |jl| will lead to p{t) 
for t 00. 

Apart from complete positivity of the time evolution, the assumptions on which the 
present work is based are the following: 

1. We work in a 2-dimensional complex Hilbert space, which we can identify with C^; 

2. we assume hermitian Lindblad operators, i.e., Aj = Aj Vj. 
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Both assumptions are crucial for the following discussions. The first one is motivated by 
the applications in particle physics, whereas the second assumption guarantees that the 
entropy S[p] = — Tr (plnp) cannot decrease as a function of time In this framework 
we will discuss the different forms of the dissipative term D[p\ used in the literature and we 
will show their equivalence. We will put emphasis on the formulation of the time evolution 
(H) in H^, where we will represent D[p\ as a positive sum over projectors onto planes. We 
will also study in detail the matrix formulation of D[p] as advocated by Benatti and 



Floreanini, e.g., in Refs. [0, |18|, and relate this formulation with the geometric version 
of D[p] as a sum of projectors. Finally, we will investigate the limit t ^ oo of the density 
matrix. 



2 Equivalent forms of the dissipative term 

We start with the original form of the dissipative term (0) and take into account that we 
confine ourselves to hermitian Lindblad operators Aj. Thus D[p] simplifies to 

Form A: D[p] = l± (A> + pA'^ - 2 A,pA,) = i ±[A„ [A,, p] ] . (9) 

The number of Lindblad terms in Eq. is denoted by r. 

Next we note that the dissipative term can be rewritten in terms of projectors Pj and 
their orthogonal complements = \ — Pj as 

Form B: D\p\ = \y^ X, (P,pP^ + P^pP,) . (10) 

The projectors Pj are non-trivial projectors in C^, which can be parameterized as 

P, = i(l + n,-a) , (11) 

where the fij are real unit vectors and a denotes the vector of Pauli matrices. The 
quantities Aj are real, positive numbers. 

Proposition 1 Forms A and B of D[p] are equivalent. 

Proof: A general hermitian 2x2 matrix Aj can be represented by 

= ^(%1 + V^^i-^) , (12) 

where aj and Aj are real numbers (Aj > 0) and fij is a unit vector. Note that the part of 
Aj proportional to the unit matrix does not contribute to D[p\, as evident from Eq. (P). 
Therefore, Aj must not be zero in order to have a non-trivial effect of Aj. Consequently, 
we are allowed to replace Aj by \J~XjPj in D[p\. Using Eq. (||) and Pj = Pj, we derive 

D[p] = lj2 A, (P,p + pP, - 2 P,pP,) , (13) 
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which can be rewritten in Form B of Eq. ([T0|) . □ 



The time evolution (|T]) is easily reformulated as a differential equation for a real 3- 
vector p by using 

1 1 . ^ ^ 

(14) 



p = ^ (1 + p ■ ct) and H = ^(^l + h- 



a 



Since for simplicity we have assumed that H is hermitian, the 3-vector h is real as well 
(for an extension to non-hermitian Hamiltonians see, e.g., Ref. |^). The new version of 
Eq. (|1]) is given by 

^ = hxp-Lp. (15) 

With Form A or B of the dissipative term and Eq. (|T^) or ([TT|), we arrive at a geometric 
version of D[p] 

Form C: Lp = - ^ Xj (p — Hj Hj ■ p) , (16) 



where the matrix L is a positive linear combination of projectors 

'P(nj) = I3 — UjU^ . 



(17) 



The projector (|T^) projects onto the plane orthogonal to rij. Clearly, Form C is equivalent 
to the the previous forms of the dissipative term. 



The dissipation matrix L in Form C can be reformulated as [|16| 



with 



Form D: L^p = 2 " ' ^Z^) 



q« G R'' (a = 1,2, 3) and A = ^ |q„|^ 

Q = l 



(19) 



Proposition 2 Forms C and D are equivalent. 

Proof: The proof of this statement amounts to a mere rewriting of the elements of L by 



j=i j=i j=i 



Then we define 



and A = A,- . 

y ^/K{nr■)a J 

The sum over the square of the lengths of the three vectors qo, is performed via 

3 3r r / 3 \ r 



0=1 a=lj=l j=l \a=l / 



EA.- 



(20) 



(21) 



:22i 



Here, we have used that the fij are 3-dimensional unit vectors. Thus we have obtained 
Form D of L from Form C. This procedure can be reversed: given three vectors Qq, G IV, 
we can use Eq. ( ^I]) to construct r unit vectors fij and positive numbers Xj. In this way, 
we gain Form C from Form D. □ 

The question arises how many terms are necessary in D[p] in the most general case. 
This question is answered by the following theorem. 

Theorem 1 The most general case is covered by three Lindblad terms in D[p], i.e., if 
D[p] is given by a sum over more than three terms, then it can be rewritten as a sum of 
at most three terms. If D[p] is formulated with the minimal number of terms, then, using 
Forms B or C, there are three distinct minimal cases referring to one, two or three linearly 
independent vectors fij with Xj > 0; we will denote these cases by an index i G {1,2,3}.[] 

Proof: Let us first assume that r > 3. Since there are only three vectors q^, we can find 
a rotation R acting on IV such that 



\ 



V 







with Qa G 



(23) 



In the most general case the set of vectors Q 
and can be parameterized by 



{Qa|a = 1,2,3} is linearly independent 



( 
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To, \ 



/ 7/ir(mi)o \ 



(24) 



The real and positive numbers pj are chosen in such a way that the vector mi, extracted 
from the first elements of the vectors Q^, is a unit vector; the same is done for the second 
and third elements. | In this case we have an index £ = 3. If Q spans a 2-dimensional 
space, we choose the rotation R such that the third elements of all Q„ are zero; if Q spans 
a 1-dimensional space, the second and third elements are taken to be zero. These two 
cases refer to index 1 = 2 and 1, respectively. Form D of the dissipative term tells us that 
D[p] is independent of any rotation R. Consequently, we arrive at 



i=i i=i 



m 



(25) 



which proves the theorem for r > 3. For r < 3 we follow analogous steps performed after 
Eq. ([23|), with the rotation R acting now on a 3 or 2-dimensional space. For r = 1 the 
procedure is trivial. □ 



■I- Obviously, in the minimal formulation of D\p\ we have r = £. 

^This procedure is analogous to the one to prove that Form C follows from Form D. 
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On the other hand, to find the minimal number of operators Aj needed for the most 
general dissipative term (0), we could start from the Gorini-Kossakowski-Sudarshan ex- 
pression (0). Transforming back to the Lindblad structure @ by using the relations 

Cki=(UcU^) and B, = J2 Ukj^j Fk , (26) 

k=l 

where c > is diagonal and U a unitary matrix, it can be seen that, in the general case, 
we obtain d"^ — 1 terms in D[p\ of Eq. (Q); if some of the elements cjj are zero, we will 
have less than d"^ — 1 terms. Applying this to c? = 2, we find at most three terms, which 
agrees with the result of the explicit calculations leading to Theorem |l|. 



3 Complete-positivity conditions on the dissipation 
matrix L 



Benatti and Floreanini fl^, |I^ parameterized the dissipation matrix L by 6 real 
constants and expressed complete positivity in the form of inequalities satisfied by these 
parameters. Thus they have the version 



(27) 





/ 


a 


b 




Form E: L = 2 




b 


a 






I 


c 


/3 


7 / 



together with 



2R 
2S 
2T 

RS>b\ 



: a + 7 — a > , 
; a + 7 — a > , 
: a + a — 7 > ; 
RT>c\ ST>(3'^] 



EST > 2bcp + Rp^ + Sc^ + Tb^ 



(28a) 

(28b) 
(28c) 



We will show now that Form E is equivalent to the forms presented in the previous section. 
First we want to put Eqs. (|28a|) , (|28b| ) , (|28c| ) into a simpler equivalent form. 



Lemma 1 Given L, there exists a symmetric matrix M such that 



1 



L = -(TrMls 



In terms of M, Eqs. 



M). 

^Sdj) , (^8d{ ) are given by 



(29) 



(i) > (« = 1, 2, 3) 

(iii) detM>0, 



(30) 



respectively. 
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Proof: For a 7^ /5 we have Mq,^ 



-2Lap- The diagonal elements of M are determined 



by 



-Laa + Lfs/j + L^^, where a 7^ /3 7^ 7 7^ a. The second part of the lemma is 



proven by plugging Eq. (^) into Eqs. ( p8a|) , ( p8b|) , (|28c|) . □ 

The next lemma will allow us to connect M possessing properties (|30D with Form D 
of the dissipative term. 

Lemma 2 A real and symmetric 3x3 matrix M has the properties of Eq. if and 
only if there exist vectors G R'^ fa = 1, 2, 3j such that 



qa ■ q/3 



(31) 



Proof: 

(<^=) This direction of the proof is quickly dealt with. Since Maa = q^ > 0, property 
(i) is valid. Furthermore, for a (3, with the Cauchy-Schwarz inequality we derive 
^ai3 = (q« ■ q/?)^ — qa q| = MaaMpi^ and property (ii) holds as well. Defining a 3 x 3 
matrix q = (qi,q2,q3), we have detM = det(g"^g) = (detg)^ > and thus property (iii). 
This completes the first half of the proof. 

(^) This direction of the proof is more involved. If all M^a were zero, then according to 
property (ii) of Eq. ( pO| ) we would have M = 0. Therefore, at least one of the diagonal 
elements of M is non-zero. Without loss of generality we assume Mn > 0. First we 
consider the case 

M11M22 -M12 > 0. 



(32) 

.2 _ .3 _ 0. 

From M12 = qi ■ q2, it follows after a sign choice that ql = M12/ \/ Mn. Taking into 
account that q2 ■ q2 = ^^22 and defining = 0, the first two vectors are given by 



Denoting the elements of by g^, we can define a vector qi by q\ = y/Mu, qf = q{ 











f 


M12/VM11 


\ 


qi = 







, q2 = 




yM22 - Myuu 






I 


) 












(33) 



With the relation of Eq. (|32|) we find that q2 is a well-defined real 3-vector. Next we use 
Mi3 = qi ■ qs and M23 = q2 • qs and obtain 



Ml 



13 



and ql 



M23 - M12M13/M: 
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(34) 



It remains to take into account q3 ■ q3 = M33. After some algebra we arrive at 

= M33 - (.3^)^ - (.3^)^ = ^^^..2 > . (35) 



MnM22 - 



The positivity follows from property (iii) of Eq. (^). Thus, gf is well-defined and we have 
proven relation (0), provided condition ( P^ holds. 
It remains to check the same for the special case 



MnM22 = M^^ , 



(36) 
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which was excluded by Eq. (|32D. From Eq. (^) we obtain M12 = r|^/M^^[M^ with 77 = ±1 
and 

' (37) 



det M = - (^^MnMss - r^y^Mia^ . 



Since det M > (see Eq. (H^)), it follows that 



M23 = r/Mi3VM22/Mn . 
With this relation and taking into account condition (PUj), it is easy to check that 



(3^ 



qi 



Mil 







(l) 




1 


M13/VM11 


\ 





, q2=V\/M22 







1 


yM33 - Mfg/Mll 




















represents a consistent choice of vectors, which fulfills equation 
proof. □ 



(39) 

). This completes the 



With Lemmata |T] and |^ we now readily see that a matrix L fulfills the conditions of 
Eqs. ( p8a| ), (|28b| ) , (|28c| ) if and only if L is given by Lap = ^{A6af3 — q^ ■ q/3), where 
q„ e R3 and A = ELi Iq^P- 

Proposition 3 Forms D and E of the dissipation matrix L are equivalent. 

Proof: According to Lemma [l|, from the matrix L we construct a matrix M with proper- 
ties (pOl). Lemma |^ tells us that such an M is represented by a matrix of scalar products 
(see Eq. (|3T|) ) and vice versa. Therefore, our statement is true. □ 



4 The asymptotic limit t ^ oo 

Now we consider the asymptotic limit of the density matrix p{t) with time evolution (|1|) 
| 20| . For this purpose we use the results of Theorem |l|, where we have also defined the 
index ^ of L. For a general study of the large time behaviour of the density matrix starting 
with expression (0) see Ref. [ p9[] . 

Theorem 2 For index £ = 2 or £ = 3, the asymptotic limit of the density matrix is given 
by 

limp(t) = il. (40) 

t— >oo ^ 

Proof: This statement is most easily proven by using Form C, Eq. (|16|), of the dissipative 
term and Eq. ([I5|) of the time evolution. Denoting by A the operator on the right-hand 
side of Eq. (^) , we have 

1 ^ 

^x = /ixx--^A,P(n,)x (41) 
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for arbitrary complex 3-vectors x. If we can show that every eigenvalue c of ^ fulfills 
Rec < 0, the theorem is proven. Let x be a normalized eigenvector with eigenvalue c. 
Then, with x^x = 1 we have 



and, consequently. 



1 ^ 

cUx = -2ih • (Rex X Imx) - - ^ AjX^P(nj)x , (42) 



1 ^ 

Re c = — AiX^P(nj)x . (43) 



Since the projectors Vijij) are positive operators and Aj > 0, we have Rec < 0. If 
Rec = 0, it is necessary that x'''P(nj)x = 1 — \nj ■ xp = Vj = 1, . . . , Therefore, the 
eigenvector x is proportional to fij Vj = This is a contradiction for £ = 2 or 3 

independent vectors fij and, indeed, every eigenvalue c has a negative non-zero real part. 
□ 

Theorem 3 For i = 1 (fii = n, Pi = P , Xi = X) we have either h \\ n or, equivalently, 
[if, P] = 0, in which case we obtain 

lim p{t) = Pp{0)P + P^p{0)P^ ; (44) 

t — ^oo 

or h ^ -n, i.e., [H,P] ^ 0, then the asymptotic limit of p(t) is the same as in Theorem^. 

Proof: We follow the same strategy as in the proof of the previous theorem. Thus either 
Rec < or the eigenvector x is proportional to n, in which case Rec = 0. Assuming 
X = r? without loss of generality, we have now ^x = h x n = en, where the eigenvalue c 
is imaginary. This equation is only soluble for h oc n, whence it follows that c = 0. The 
relation /i oc ri is equivalent to [H, P] = 0. In this case we can write the Hamiltonian as 
H = HP + h'P-^. Decomposing an arbitrary density matrix p as 

p = po + Pi with po = PpP + P^pP^ and pi = PpP^ + P^pP , (45) 

we find 

^ = and ^ = [~i{h - h') - A/2] PpP^ + [i{h - h') - A/2] P^pP . (46) 

Therefore, we arrive at 

po(t) = po(0) and pi(t) = e^-'if^-'^'^-mt Pp(o)P^ + e^iih-h'yx/2]t p±p(o)P . (47) 

Since lim(_»oo Pi{t) = 0, the theorem is proven. □ 

The different limits of p(t) discussed in Theorems ^ and ^ have been noticed in 



Refs. H, ||. In the case of £ = 1 and [H,P] = 0, the limit (||) of the density 
matrix has the form p = pP + (1 — p)P'^ with < p < 1; all density matrices obeying 
dp/dt = have this form for [H, P] = 0. For ^ = 1 and [H, P] ^ 0, and for ^ = 2, 3, the 
unique density matrix which is time-independent is proportional to the unit matrix. 
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5 Summary and Conclusions 



In this paper we have considered the quantum-mechanical time evolution of a 2 x 2 
density matrix p{t), where the von Neumann equation is modified by a dissipative term 
D[p], which, therefore, must be of the Lindblad type or, equivalently, of the Gorini- 
Kossakowski-Sudarshan type in order to assure a completely positive time evolution. We 
have, furthermore, assumed that the Lindblad operators Aj are all hermitian which ensues 
that the entropy is non-decreasing with time. Our starting point conforms with many 
applications of the open-systems approach in particle physics. 

We have discussed five equivalent forms of D[p\ which all have their merits depending 
on the problem considered. We have put particular emphasis on the time evolution in the 
form of Eq. (|T^), where the density matrix and the Hamiltonian are represented by real 
3- vectors, and Form C, Eq. (|16]), of the dissipative term, where the dissipative term is a 
positive linear combination of projectors onto 2-dimensional planes in R^. 

We have studied the question of the minimal number i of Lindblad terms needed in 
order to reproduce a given D[p\ and formulated the result in Theorem ^ the proof of 
this theorem represents at the same time a procedure how to determine i in practice. 
An other procedure would be to start with the Gorini-Kossakowski-Sudarshan expression 
from which it also can be seen that D[p] can be generally decomposed into 3 terms. 

We have also connected the approach where the dissipative term is given by a matrix 
L specified by the conditions ( p8a| ), ( p8b| ), ( |28c| ) with the geometric picture of D[p] as 
given by Form C, Eq. ([T6|). Again, the proof which shows the equivalence between the two 
approaches, given by Lemmata || and ||, indicates a practical way to obtain the projectors 
Vifij) ([T7|) associated with D[p\. 

Finally, we have presented a general discussion of the limit t ^ oo of p(t), where the 
usefulness of Form C of D[p\ was exemplified. 
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